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Abstract 

We demonstrate that delta-doping can be used to create a dimensionally confined region of 
metallic ferromagnetism in an antiferromagnetic (AF) manganite host, without introducing any 
explicit disorder due to dopants or frustration of spins. Delta-doped carriers are inserted into a 
manganite superlattice (SL) by a digital-synthesis technique. Theoretical consideration of these 
additional carriers show that they cause a local enhancement of ferromagnetic (F) double-exchange 
with respect to AF super exchange, resulting in local canting of the AF spins. This leads to 
a highly modulated magnetization, as measured by polarized neutron reflectometry. The spatial 
modulation of the canting is related to the spreading of charge from the doped layer, and establishes 
a fundamental length scale for charge transfer, transformation of orbital occupancy and magnetic 
order in these manganites. Furthermore, we confirm the existence of the canted, AF state as was 
predicted by de Gennes [P.-G. de Gennes, Phys. Rev. 118 141 1960], but had remained elusive. 
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In semiconductors, delta-doping [l| has led to record high mobilities in two-dimensional 
electron gases, spawning a number of fundamental discoveries and many novel applications. 
In the complex oxides, similar doping strategies may in principle be used to create two- 
dimensional analogs of the many collective phases found in these materials 

as. 

In this 

work, using AF Lai-^Sr^MnOs as the host material, we have devised a delta-doping strategy 
for locally tailoring the strength of the F double-exchange interactions relative to the AF 
superexchange, creating a quasi two-dimensional region of enhanced moment. We do this 
without introducing a random disorder potential due to dopants, or the frustration inher- 
ent to some AF/F interfaces. Our approach is unique to systems where double-exchange 
(DE) and superexchange (SE) interactions compete, and cannot be realized in conventional 
semiconductors or metals. 

In the seminal work of de Gennes [4 1 , it was recognized that when the insulating AF parent 
compound of a manganite is doped with carriers, competition between DE and SE causes 
the AF spins to cant, in proportion to the doping. In real materials, additional effects such 
as the Jahn- Teller (JT) interactions, on-site Coulomb repulsions (Mott-Hubbard effects) 
and orbital ordering instabilities, all favor localizing the carriers and compete with DE. As 
a result, the transition between an AF insulator such as LaMnO"3 to a F metal such as 
Lao.8oSro.2oMnO"3 does not follow de Gennes' phase diagram but rather takes place through 
a mixed-phase region for intermediate values of x pj. A tendency toward phase separation 
is believed to be a fundamental property of the double-exchange manganites 6] and may 
thwart the appearance of a homogeneous, canted AF phase Q. 

The situation is quite different if one starts with Lai_ x Sr x Mn03 for x > 0.5, which 
is also AF. The Mn e g electrons can go into one of two partially filled bands made of 
degenerate d 3z 2_ r 2 and d x 2_ y 2 orbitals. For 0.5 < x < 0.7 this degeneracy is removed by the 
formation of a highly anisotropic A-type AF order (planes of F spins that are mutually AF), 
with the electrons mainly occupying the d x 2_ y 2 orbital states 8|. The transport is nearly 
metallic in-plane due to DE. However, the out-of-plane resistivity is orders of magnitude 
higher since the d 3z 2_ r 2 states are unoccupied and SE dominates in this direction |Sj. Upon 
doping with electrons toward x < 0.5, the d 3z 2_ r 2 orbitals begin to be occupied, DE begins 
to act in the out-of-plane direction, and the material transforms into a 3-dimensional F 
metal. We show in this work that by delta-doping an AF manganite SL with electrons, 
we locally enhance the DE which cants the spins away from AF alignment, resulting in a 
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highly modulated magnetization. Our theoretical investigation shows that a JT distortion 
is instrumental in stabilizing this canted spin structure lying in the midst of the crossover 
from a 2D antiferromagnet to a 3D ferromagnet near x — 0.5. 

We synthesized cation-ordered analogs of Lao.sSro.sMnOs by alternating single unit cell 
layers of LaMnOs (LMO) and SrMnC>3 (SMO) with atomic layer precision onto SrTi03 
substrates using ozone-assisted molecular beam epitaxy. Neutron diffraction measurements 
on these (LaMn03)i/(SrMn03)i superlattices confirmed the A-type AF phase, and mag- 
netometry measurements using a superconducting quantum interference device (SQUID) 
showed a net F moment near zero at low field n]]. Here, we present delta-doped SLs with 
nominal compositions x = 0.44 and 0.47. The x = 0.44 SL was made by alternating unit 
cell layers of LaMnC>3 and SrMn03 and inserting an additional LaMn03 layer for every 
four LaMnOs/SrMnOs bilayers ll). This sequence was repeated 9 times (9 supercells), 
to form the superlattice [(SM01/LM01)x4,LM01]x9. A La .5 6 Sr .44MnO 3 alloy film with 
a thickness of 81 uc was made for direct comparison. The x = 0.47 SL was made by in- 
serting an extra LaMnOs layer after every 9 LaMnOs/SrMnOs bilayers, repeated 4 times: 
[(SM01/LM01)x4,LMOUSM01/LM01)x5]x4. Structural characterization of the SLs is 



discussed elsewhere 
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111 ]. By periodically inserting the extra LaMnC>3 unit cell layer, we 



have delta-doped the MnC>2 planes in the vicinity of the inserted LaMnC>3 unit cell with 
electrons. Our investigation of these SLs by neutron scattering techniques was aimed at 
determining the influence of these delta-doped charges on the A-type spin structure. 

We performed neutron diffraction measurements to probe the spin structure in the x = 
0.44 SL and alloy. The A-type order is manifested by a structurally-forbidden (0 |) Bragg 
peak, shown in Fig. [TJi, signifying an AF spin alignment with a periodicity of 2c normal to 
the film plane (c is the out-of-plane lattice parameter). From the full width at half-maximum 
of the peak, a magnetic coherence length of ~28.5 nm was determined, nearly the entire film 
thickness (30.9 nm). The coexistence of both a (0 |) neutron diffraction peak and a net F 
moment from M(H) data ll| suggests a canted AF structure, for which the spins within a 
single plane are ferromagnet ically aligned and the spins on adjacent planes are canted at a 
relative angle 8 Z < 180°, where z /2 is the angle between each spin and the applied magnetic 
field direction. The temperature (T) dependence of the peak intensity in Fig. [lb indicates 
that the A-type order in the SL disappears at ~260 K. 

To determine the depth-dependent magnetization profile, we performed polarized neutron 
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FIG. 1: (a) (0 ^) neutron diffraction peak for the x = 0.44 SL. The lines are a guide to the eye. 
(b) T-dependence of the integrated intensity of the diffraction peak for the x = 0.44 SL and alloy, 
measured while increasing T in zero field, after cooling in zero field. 

reflectometry (PNR) measurements, which detect variations of in-plane magnetization as a 
function of sample depth along the surface normal. The non-spin-flip reflectivities, R ++ and 
R shown in Fig. [2^ for the x = 0.44 SL, are sensitive to the component of magnetization 
parallel to the applied field direction. Because the nuclear scattering length densities of SMO 
and LMO are nearly identical, only magnetic variations in the layer profile are detectable. 
The key feature of interest here is the Bragg peak at q = 0.18 A -1 , indicating that the 
magnetization is in fact modulated with a periodicity that matches the supercell periodicity 
of the SL. 



To fit the PNR data using the co_refine routine 12] . a model consisting of 9 supercells was 
used, where each supercell consisted of a high moment sublayer and a low moment sublayer. 
The sublayer thicknesses, magnetic roughnesses and magnetic scattering length densities (a 
quantity directly proportional to the magnetization) of the seven interior supercells were 
fit while constrained to be identical, whereas the topmost and bottommost supercells could 
vary independently. For this data and model, there are several sets of possible fit parameters 
that yield the same lowest x 2 value to within 1%. Two such profiles are shown in Fig. [2b. 
Even though the highest (lowest) moment in the modulation period, M max (M min ) can each 
vary by as much as 0.25 hb, we can conclude with certainty across all sets of fit parameters 
that 1) the period of the magnetic modulation matches the structural supercell to within 
0.1 nm; 2) the high moment sublayer contains the extra LMO layer and 3) the magnitude 
of the modulation is quite large, with M max /M m i n = 3 ± 0.3. 

The x = 0.47 SL has a thicker supercell with one extra LMO layer for every 19 uc, such 
that two Bragg peaks are detectable, as shown in the PNR spectra in Fig. [2b. With this 
additional information on the magnetic modulation compared to the x = 0.44 SL data set, 
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FIG. 2: (a) The non-spin-flip channels (R ++ and R~~) of the PNR for the x = 0.44 SL at 
T = 120 K and with H = 675 mT applied in the plane of the film, has a Bragg peak signifying that 
the magnetization is modulated. The solid lines are the fit to the data. The inset shows the PNR 
of the x = 0.44 alloy film having no Bragg peak, (b) Two possible magnetic profiles as a function 
of film depth from the fit to the PNR spectra of the x = 0.44 SL. The location of the LMO (pink) 
and SMO (green) layers in the SL and the SrTiO"3 substrate (gray) are shown, (c) PNR for the 
x = 0.47 SL with H = 815 mT and T = 120 K. The Bragg peaks are more apparent in the inset, 
which shows the same data multiplied by g 4 to compensate for Fresnel decay, (d) The magnetic 
depth profile of the x = 0.47 SL from the best fit. 
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the fitting routine produces a single set of best fit parameters at the lowest x 2 f° r the x = 0.47 
SL. A model similar to the one previously described, containing alternating high moment 
and low moment sublayers, was used to fit this data. The moment is again highly modulated, 
alternating between a high moment region across 6 uc (the full-width half-maximum of the 
peaks in Fig. [2U) containing the extra LMO unit cell and having M max = 2.2/ig, and a 
low moment region spanning 13 uc having M m j n = 0.7^b- The gradual transition between 
M max and M min occurs over 4 uc, and as we shall show, this occurs by varying the canting 
angle. Using M to t = (4 — x)hb as the magnitude of the moment of a Mn 3+//4+ spin and 
M = M tot cos(8 z /2), the measured modulation of M corresponds to a modulation of canting 
angle between 9 Z = 103° and 157°. These fit parameters also work for the x = 0.44 SL, where 
the high moment region extends across 6 uc with M max = 2.2/ig, and the low moment region 
extends across 3 uc with M m i n = 0.7/zb, as shown in Fig. [2b (black line). In comparison, the 
PNR measurement of the x = 0.44 alloy film (inset of Fig. [2k.) shows a net moment without 
any modulation (no Bragg peak). Thus, the spin canting in the alloy is uniform, and the 
modulated magnetization in the x = 0.44 and 0.47 SLs occurs due to the delta-doping. 

To investigate this canted spin structure further, we carried out diffraction measurements 
with polarized neutrons and polarization analysis in the vicinity of the AF peak for the 
x = 0.44 SL, as shown in Fig. [3k. The peak in spin-flip intensities, I + ~ and at 
q = 27r/(2c) = 0.82 A -1 results from the canted AF spins with each sublattice having a 
component of magnetization pointing perpendicular to the applied field direction, with a 
period of 2c. Moreover, the superlattice peak at q — 0.64 A -1 signifies that the canted spin 
structure is modulated with a period of 9 uc, consistent with the supercell period and the 
periodicity of the F component revealed by the low q PNR measurements [l3]. These peaks 
in J H and I h directly confirm the canted, modulated spin structure. Fig. [3b shows the 
canted, modulated spin structure as deduced by the PNR and diffraction measurements. 

In order to explore theoretically this canted spin state, we consider the standard two- 
orbital model for manganites [7] with the Hamiltonian H = H^e + H$e + Hjt U|. The 
double-exchange term H^e, which arises from a combination of kinetic energy and strong 
Hund's rule coupling, favors the F state. The superexchange term H$e is the Heisenberg 
interaction between t 2 g-spins favoring AF coupling. The Jahn- Teller term Hjt represents 
the coupling of the e g electrons to the lattice distortions. 

Before understanding the canted magnetic states, it is essential to understand the stability 
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FIG. 3: (a) Diffraction with polarized neutrons and polarization analysis, measured at T = 120 K 
in a 820 mT field. Peaks in the spin-flip intensities at high q directly confirm the canted, modulated 
spin structure. The inset shows the non-spin flip intensities, (b) Schematic of the canted, modulated 
spin structure, showing 4 supercells (9 uc each) of the x = 0.44 SL. Pink and green arrows represent 
the F spin alignment within the MnO"2 planes of the LMO and SMO layers, respectively. The dotted 
blue (purple) lines are the projection of the spins onto the yz-plane (xz-plane) and signify the 
magnetic potential measured by neutron diffraction (PNR). The dotted red lines are the projection 
of the spins onto the xy-plane, showing the canting angles. The black lines aid the eye in identifying 
the spin modulation. 

of A-type AF order near x = 0.5 in manganites. It is an outcome of a competition between 
the ferromagnetic H DE = J F cos{9ij/2) and the antiferromagnetic H SE = JsJ2 C0S (^ij) 
terms, where 9ij denotes the relative angle between spins Sj and Sj, and the summation is 
over all nearest neighbors. A-type order requires the F term to dominate in the xy-plane 
and the AF term to be stronger along the z-axis. The presence of a uniform JT distortion 
favors the occupation of d x 2_ y 2 orbitals, which leads to a strong F interaction in the xy- 
plane and a much weaker one along the z-axis, as seen in Fig. [4^ (Jp ~ p z t z , where p z is 
the d 3z 2_ r 2 orbital occupancy and t z is the hopping integral in the z direction). This leads 
to A-type AF order for x > 0.5 when the SE coupling J$ >> Jp. Upon reducing x so that 
filling of the e g band continues, the doped electrons begin to fill the higher e g states (the 
ds z 2_ r 2 states), so that Jp increases and competes with J5 while the in-plane order remains 
unaffected, as shown in Fig. Hk- Such a competition between F and AF interactions along 
the z-axis can lead to canted spin states. From a simple energy minimization dE/d6 z = 
for E = —J F cos(9 z /2) + J s cos(9 z ), one finds that the canting angle 9 Z = 2 arccos( J F /4J S ) 
if Jp < 4J5, and 9 Z = (F alignment) when Jp > 4J^. Note that within this scenario one 
needs Js = 00 to obtain 9 Z = 180° (AF state). In reality there is a feedback effect of the 
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FIG. 4: (a) The effective F exchange constant arising from the DE interaction, shown here for a 
fixed JT distortion and a ferromagnetic spin state for all x. The presence of a JT distortion leads 
to a difference in orbital occupancy, which in turn leads to anisotropic exchange parameters. A 
typical value of the AF coupling, J$ = 0.14t is shown as a dashed line, (b) Doping dependence 
of the canting angle Z for various values of J5. The canting angle is obtained by minimizing the 
total energy at a given x, where both the JT distortion and 9 Z are allowed to vary. 

spin states on the orbital occupancies that causes the variation of 9 Z to be more abrupt. 
When the spins cant away from F alignment, this causes occupancy of the d^ z 2_ T 2 orbitals 
to decrease, which in turn causes Jp to decrease further, making the spins cant even more. 
The cost in elastic energy in creating the lattice distortion with changes in orbital occupancy 
(c-axis shrinks with lower p z ) prevents the transition from being a step function. 

In order to test this explicitly, we minimize the total energy of the Hamiltonian for 
variational parameters 9 Z and JT distortions, hence allowing for this feedback effect. At a 
fixed x, JT distortions lead to a reduction in electronic energy but there is also an elastic 
cost associated with making these distortions [ljj. The optimum value of the distortion is 
selected from the above competition. This generates an x dependence on the magnitude of 
the distortions and therefore an additional x dependence on the values of Jp. The resulting 
optimum angle 9 Z is plotted in Fig. H)d for a few values of J5. The sharp change in 9 Z from 
180° to 0° for Js = 0.10, 0.11 in the calculation does not allow for layers at intermediate 
x with 0° < 9 Z < 180°. This implies that it is not possible to stabilize a finite canting 
angle for these parameters. However, for J$ > 0.12 and for x < 0.5, the variation in 9 Z 
with x is more gradual and layers with a finite canting angle are stable. Therefore, a small 
variation in the doping level can lead to significant changes in canting angle and thus the 
net F moment. The large modulation in moment as a function of depth in these delta-doped 
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SLs is a realization of this scenario. Attempts to realize the canted AF state near x = as 
proposed by de Gennes and near x = 0.5 have typically been obscured by phase segregation 
into F and AF domains [jjj], which is mitigated in our digital synthesis approach of charge 
doping without disorder. We note that a modulation in orbital occupancy is expected to 
accompany this modulated spin state [151 ]. 

In conclusion, we have employed delta-doping to tailor the magnetic structure following 
the ideas of de Gennes. Localized regions of enhanced F exchange were created in an 
otherwise AF structure, where the transition between the high moment region and the low 
moment region occurs via a continuous variation of canting angle that is governed by the 
spreading of charge near the delta-doped layer. Furthermore, the PNR technique enables us 
to determine a length scale of 6 uc for the region of enhanced moment that is induced by the 
delta-doped layer. The result that the SL magnetization is highly modulated is evidence that 
the doped charges are not completely delocalized over the entire SL. Charge spreading out 
to 3 uc from both sides of the delta-doped layer is consistent with previous theoretical and 
experimental works that had inferred a length scale for charge spreading via more indirect 



means 
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171 ] . We believe this to signify a fundamental length scale for the spreading of 



charge normal to the layers. 
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